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We show that the energy threshold for nuclear recoils in the XENON 10 dark matter search data 
can be lowered to ~ 1 keV, by using only the ionization signal. In other words, we make no 
requirement that a valid event contain a primary scintillation signal. We therefore relinquish 
incident particle type discrimination, which is based on the ratio of ionization to scintillation 
in liquid xenon. This method compromises the detector's ability to precisely determine the z 
coordinate of a particle interaction. However, we show for the first time that it is possible to 
discriminate bulk events from surface events based solely on the ionization signal. 
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1. Introduction 

The XENON 10 detector [1] is a liquid xenon time-projection chamber with a an active tar- 
get mass of 13.7 kg (15 cm height and 10 cm radius). It was designed to directly detect galactic 
dark matter particles which scatter off xenon nuclei. Typical velocities of halo-bound dark matter 
particles are of order 10~ 3 c. This leads to the prediction of featureless exponential recoil energy 
spectra for spin- independent elastic scattering of dark matter particles on a xenon target [2]. Typ- 
ical nuclear recoil energies are few keV. A particle interaction in liquid xenon results in a prompt 
scintillation signal (SI) and an ionization signal (S2). The generation and collection of these sig- 
nals are discussed in detail in [1]. A key feature of the XENON10 detector is the event-by-event 
discrimination between incident particle type, based on the ratio S2/S1. 

The energy threshold of previously reported XENON10 data [3-5] depends on the primary 
scintillation efficiency of liquid xenon for nuclear recoils ££ e ff [6, 7], and for a conservative as- 
sumption of the energy dependence of =S? e // [6], is about 5 keV. This energy threshold is dictated 
by the collection of primary scintillation photons following a nuclear recoil in the xenon target. But 
it is possible to obtain a lower energy threshold from the existing XENON 10 data, using only the 
ionization signal. We discuss an analysis of XENON10 dark matter search data, with the energy 
scale set by the detected ionization signal and no requirement that valid events contain both an 
SI and an S2 signal. This requires a compromise on two important aspects of the detector per- 
formance: the ability to precisely reconstruct the z coordinate of a particle interaction (which is 
normally obtained by the time delay between SI and S2 signals), and the discrimination between 
incident particle types. 

Calibration of the ionization energy scale for nuclear recoils in liquid xenon is a pre -requisite. 
Direct measurements exist in the literature [6, 8] (and are shown in Fig. 2b), but do not extend to 
the smallest signals observed in our detector. We describe a new measurement of the ionization 
yield of liquid xenon for nuclear recoils, down to a nuclear recoil energy of about 1 keVr. We then 
describe a method for making an approximate determination of the z coordinate of the interaction. 
The method exploits the measured width of the S2 signal, which is broadened by electron diffusion 
as a cloud of electrons drifts across the liquid xenon target. The precision is low (~cm), but is 
sufficient to discriminate bulk events from surface events. This is important because the latter are 
more likely to result from radioactive background. 

2. Ionization yield of liquid xenon for nuclear recoils 

The ionization yield J2 y was obtained directly from in-situ neutron calibration data. The neu- 
tron calibration experiment and general analysis technique are both described in a previous study 
by the XENON 10 collaboration [9]. Note that here we circumvent the. false single scatter pathol- 
ogy described in [3], since our calibration technique relies only on the S2 signal. The S2 spectrum 
for single elastic neutron scatters is shown in Fig. la, with la uncertainty. The S2 signal was 
scaled to an absolute number of electrons via the measured 24 photoelectrons per liquid electron 
[1]. Only the most basic data quality cuts were applied, along with a fiducial cut r < 8 cm. No 
z fiducial cut was applied because a substantial fraction of events at low energy have no S 1 (and 
hence have indeterminate z). The fiducial target volume used in this calibration thus differs from 
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electrons extracted from liquid xenon E nt [keVr] 

Figure 1: (a) Spectrum of the number of electrons extracted from single scatter elastic nuclear recoil inter- 
actions in the fiducial target defined by r < 8 cm, with 1 a statistical uncertainty. The best-fit Monte Carlo 
is shown as a continuous stair-step. The energy values of the 9 spline points are identical in both plots. 
Shown in light gray is the spectrum of events in which an SI signal was found in the 80 jJ.s before the S2. 
Approximate keVr-equivalent is indicated along the top. (b) The ionization yield from nuclear recoils £} y 
from the present work (#), with la statistical uncertainty. Systematic uncertainty is shown along the axis. 
Also shown are data from [6] (*), [9] (V) and [8] (O and ■ ). 

previous analyses [4, 5, 3, 9] in that it considers the full Az = 15 cm active target (8.6 kg target 
mass), rather than just the central Az = 9.3 cm (5.4 kg target mass). 

With no self-shielding liquid xenon above or below the chosen xenon target volume, it is 
important to take account of possible gamma and beta background contamination in the nuclear 
recoil data. This background is assessed by counting the number of events within ±2a of the 
electron recoil centroid, in the region around 30 — 35 keVee. The chosen energy window is above 
the tail of the elastic nuclear recoil distribution, and below the 40 keVee inelastic scatters [9]. The 
electromagnetic background spectrum below about 50 keVee is predicted and observed to be flat 
with energy. Therefore, we were able to subtract the resulting background prediction of about 0.75 
counts/electron. Considering the number of nuclear recoil events in the data sample shown in Fig. 
la, this is a very small correction. 

In comparing the measured energy spectrum with the Monte Carlo predicted spectrum, the 
energy resolution in the electron signal was assumed to be Poisson in the number of detected 
electrons; account was also taken of the 20% la width of the single electron distribution [10]. The 
S2-sensitive trigger threshold for the neutron calibration data had full efficiency for events with at 
least 182 photoelectrons [1], or 7.6 S2 electrons. This is valid for r < 3 cm, and there is a slight 
radial dependence to the trigger such that by r = 8 cm the efficiency is unity for events with greater 
than 11 S2 electrons. As a result, events with S2 < 12 electrons were not used in the calibration 
analysis. No attempt was made to model the S2 trigger efficiency, as is apparent in Fig. la from 
the discrepancy between data and simulation below ~ 8 electrons. 

The ionization yield J2 y was modeled as a continuous cubic Hermite spline interpolation as 
in [9], and a maximum likelihood comparison was applied to find £} y (E nr ) as shown in Fig. lb 
(#, with la statistical uncertainty). No constraints were applied to the spline other than the re- 
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quirement that £! y = at E„ r = 0.1 keVr. The exact minimum recoil energy that can result in 
a single detectable electron is not known, but is limited by the ionization potential of xenon and 
by nuclear recoil quenching [11]. Our chosen boundary condition is conservative, and our results 
are insensitive to a x5 shift in this value. The spline points were chosen at fixed values of keV 
nuclear recoil energy (keVr), and the results do not depend on the location of the spline points (as 
long as there are ~ 10 points spanning the full energy range). The best fit shown in Fig. la has 
X 2 /d.o.f. = 414/425. For reference, the ionization yield of 122 keV gammas from a 57 Co source 
were found to have an ionization yield of 46.5 ±7.0 (stat.) electrons/keV, for events with radial 
positions in the range 8.5 < r < 9.0 cm. 

Two sources of systematic uncertainty are indicated as vertical bars for each spline point, along 
the *-axis. These arise from an assumed ±10% uncertainty in the single electron calibration (left 
bar), and from uncertainty below 2 MeV in the spectrum of initial neutron energies E„ from the 
AmBe neutron source [12] (right bar). The neutron energy spectrum was taken from Fig. 5 of 
[12], and the uncertainty was parameterized as 1 ±exp(— E„ — g) for E„ < 2 MeV. The effect of 
this conservative assumption on the Monte Carlo nuclear recoil energy spectrum was a change of 
about ±15% in bin counts at 1 keVr. A third source of uncertainty, arising from the Xe(n,n)Xe 
elastic cross-section data [13], would appear almost point-like in Fig. lb and is not shown. Our 
£ly results were found to be essentially unchanged if the assumed energy resolution was varied by 
±25%, and remained very similar if the fit range was instead truncated at either 10 or 15 electrons. 
The agreement between this measurement and that of [6] (Fig. lb, *) is quite good above 6 keVr. 
A possible reason for the rapid rise in J2 y below 6 keVr obtained in [6] is discussed in [14]. The 
previous work [9] (Fig. lb, V) shows systematically higher values of £2 y below about 25 keVr, 
possibly due to the systematic effect of false single scatters mentioned at the beginning of Sec. 2. 
It is also worth emphasizing that the previous work [9] relied on the SI signal, and therefore had 
very limited sensitivity to recoil energies smaller than about 5 keVr. 

We have argued that our nuclear recoil data sample does not contain a significant electromag- 
netic background contamination, and that what little contamination exists can be safely subtracted. 
This argument is weakened for events in which no S 1 signal was detected, since they have an in- 
determinate z coordinate for the scatter vertex, and also an indeterminate discrimination parameter 
(S2/S1). It is clear from Fig. la that most events with no SI have S2 < 40 electrons. A compelling 
piece of evidence that such events are in fact elastic nuclear recoils is that their spectral shape is 
similar to that obtained for elastic nuclear recoils with known SI. In contrast, the population of 
events with no SI in the gamma calibration data shows the expected flat S2 spectral shape in the 
range S2 < 40 electrons. However, this does not fully exclude the possibility of an additional low- 
energy background. The most likely origin of such a background, if it exists, would be low-energy 
gamma or beta scattering at the liquid surface. Bulk events would be excluded by virtue of the 
~/im range of ~keV particles in liquid xenon. 

About 95% of events at the liquid surface are measured to have 0.13 < Osi % 0.23 fis, as can 
be seen in Fig. 2a. The distribution of a si for bulk events, and events with indeterminate drift 
time are both roughly Gaussian. The former shows a 14.7% ±2.2% excess of events in the region 
0.13 < Os2 < 0.23 /is. This places an upper limit on the fraction of events with no S 1 which might 
be unaccounted surface background events, rather than bulk elastic nuclear recoils. The effect of 
such a background (assumed flat in the range 7 — 40 electrons) on our energy calibration would be 
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Figure 2: The S2 pulse width O52 obtained by Gaussian fit, for single scatter events with 7 < S2 < 100 elec- 
trons in the region r < 3 cm. The liquid-gas interface is at z = cm. Events with indeterminate z coordinate 
due to absent SI are shown at z=15 cm. The band mean is indicated by O, with Id width; horizontal bars 
indicate bin width. A 0.4 fis section of two typical S2 ~ 12 electron events are shown inset (top). In the 
event at z = 0.3 cm, the SI is visible 1.6 jis before the S2; in the event at z — 13.7 cm, the SI is far off scale. 

to increase J2 y at 8 keVr and 16 keVr by about 0.2 electrons/keVr, and by half that amount at other 
energies E nr < 32 keVr. This is a very small effect and is not indicated with the other uncertainties 
in Fig. lb. 

3. Obtaining the z coordinate of a particle interaction from the S2 width 

A cloud of ionized electrons resulting from a particle interaction drifts through the liquid xenon 
target at about 0.20 cm /is -1 [15]. As it drifts, its spatial extent broadens due to diffusion. The 
amount of diffusion broadening is reflected in the width of the S2 pulse. In this way, the width of 
S2 signals depends weakly on the z coordinate of the scatter vertex. The width Os2 was obtained by 
Gaussian fit to each S2 pulse, and is shown in Fig. 2a for events with r < 3 cm. Events with no S 1 
pulse do not have a precisely defined drift time, and are shown (arbitrarily) at z = 15 cm. A 0.4 /is 
section of two typical events with S2 ~ 12 electrons are shown inset in Fig. 2a, corresponding to 
nuclear recoils at the extrema (in z) of the active xenon target. The height of the S2 pulse in the 
z = 0.3 cm event is about 10 mV referred to the amplifier input. The z dependence of Gs2 is clearly 
weak, however it is sufficient for us to discriminate events which occur at the surface (z = cm and 
z = 15 cm) from those which occur in the bulk of the liquid xenon target. 

4. Dark matter search data 

Figure 2b shows the differential event rate from a 12.5 live day exposure of the XENON 10 
detector, obtained between August 23 and September 14, 2006. This data set is distinct from the 
previously reported [3, 4] dark matter search data. The most notable difference is that during this 
time, the detector was operated with an S2-sensitive trigger threshold at the level of a single ionized 
electron. Event selection criteria were applied as follows. Valid single scatter events were required 
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to have only a single S2 pulse in the event record. The efficiency for making this selection is very 
high, considering the robustness of the S2 signal (as shown in Fig. 2a). A fiducial cut defined 
by r < 3 cm was imposed, giving a target mass of 1.2 kg. This central region features optimal 
self- shielding by the surrounding xenon target. No explicit z cut was made since events were not 
required to contain an S 1 signal. Events in which an S 1 signal was found were required to have 
logio(S2/Sl) in the ±3a band for elastic single scatter nuclear recoils [3]. Events in which no SI 
signal was found were assumed to be low-energy nuclear recoils and were retained. The fraction of 
events with no S 1 is given by the ratio of the histograms in Fig. 1 for the nuclear recoil calibration 
data. For the dark matter search data, the fraction is approximately x3 lower. This behavior is 
expected for electron recoils, based on the observed ratio of logio(S2/S 1) for electron and nuclear 
recoils [3]. The combined acceptance of the two cuts is e > 0.99, for events in the 1.2 kg central 
region. 

We then used the S2 width Osi to discriminate bulk events from edge events, as described in 
Sec. 3. In order to ensure complete rejection of events at the liquid xenon surface (at z = cm), we 
set the lower bound of the cut at Osi > 0.23 /xs. The upper bound of the cut was set at Cs2 < 0.30 jJ.s. 
Considering the total (~Gaussian) distribution of recorded Os2 values, these bounds correspond 
approximately to the region between /I and /I + 2a. It is clear from Fig. 2a that this only partially 
targets events from the bottom of the active region of the detector (near z = 15 cm). This choice 
does not compromise surface event rejection, considering that below the cathode grid which defines 
z = 15 cm are an additional 1.3 cm of self-shielding liquid xenon. The acceptance of the S2 
width cut for single scatter nuclear recoils is mildly energy-dependent, rising monotonically from 
£ = 0.39 for events with S2 = 7 electrons, to e = 0.44 for events with S2 > 40 electrons. 

The S2 spectrum of all single-scatter events that passed these cuts is shown in Fig. 2b (O). 
Vertical bars indicate statistical uncertainty, and horizontal bars indicate bin width. The count 
rate is adjusted for the total acceptance fraction e. The lowest energy events remaining above an 
analysis threshold of 7 electrons have S2 signals of 8, 17, 18 and 29 electrons. An additional 4 
events were found with 30 < S2 < 44 electrons. For a higher statistics comparison, the spectrum 
of all single-scatter events within an 8 cm radius (target mass 8.6 kg) is also shown (□). No S2 
width cut was applied to the data in that case, and the statistical uncertainty is smaller than the 
data points. Both spectra are essentially flat above S2 = 7 electrons, as would be expected from 
a Compton scatter background. Note that in the calibration analysis of Sec. 2, a lower bound of 
7 electrons in the S2 signal was motivated by the trigger efficiency during the neutron calibration. 
In the dark matter search data the trigger threshold is at the level of a single electron in the S2 
signal. However, we retain the 7 electron lower bound in analyzing this data because (i) we can 
most accurately determine acceptance for nuclear recoils above this value, and (ii) the electronic 
noise increases significantly below S2 ~ 6 electrons, as shown in Fig. 2b. 

5. Summary 

The measured number of S2 electrons can be scaled to nuclear recoil equivalent energy via 
the £2 y curve shown in Fig. lb. A conservative choice would be the —la contour, which is ap- 
proximately £} y k, 4 electrons/keVr below E m - = 20 keVr. As mentioned in Sec. 2, the detector 
resolution for S2 signals depends primarily on Poisson fluctuations in the number of detected elec- 
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trons, with an additional component due to instrumental fluctuations. This is discussed in detail in 
[14], and for higher energy signals in [1]. So as not to overstate the energy resolution, we suggest 
a parameterization ^(keV) = 1 /2vkeV which follows the Poisson component only. This is valid 
for the — la contour of J2 y . We have shown that it is possible to give up the usual incident particle 
type discrimination based on logio(S2/Sl), and analyze the dark matter sensitivity of XENON10 
using only the S2 signal. The advantage of this analysis appears to be an increased sensitivity to 
light (< 10 GeV) dark matter particles, due to the significantly lower energy threshold. For larger 
particle masses, the usual analyses [3, 4] offer superior sensitivity. Dark matter exclusion limits 
obtained from the present work should offer substaintial constraints on recent interpretations [16- 
21] of the excess low-energy events observed by CoGeNT [22] and CRESST-II [23], as well as the 
DAMA modulation signal [24]. 
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